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O | Abstract 

We construct an electromagnetic field whose scalar potential is a pulsed- 
beam wavelet ^> (an analytic continuation of a classical Huygens wavelet). 

S The vector potential A is determined up to three complex parameters by 
requiring that (a) it satisfies the Lorenz gauge condition with (b) its 
current density is supported on the same disk T> as the charge density of 
$ , (c) it is axisymmetric, and (d) it has the same retarded time depen- 
dence as By choosing one of the parameters in A appropriately, the 
00 electromagnetic field generated by the four-potential (A, can be made 

C*"") mill, meaning that E 2 = B 2 and E ■ B = 0. We call such fields coherent 

because upon being radiated, they do not loiter around the source, gener- 
ating electromagnetic inertia, but immediately propagate out at the speed 
of light. The coherent EM wavelets define a twisting null congruence of 
light rays in Minkowski space, which we show to be identical to the Kerr 
y—i congruence associated with the Kerr-Newman metric. The latter repre- 

sents a black hole due to a time-independent charge-current density on a 
massive disk T> spinning at the angular velocity c/a, where a is the ra- 
• i— I dius of T>. By contrast, our coherent wavelets are electromagnetic pulsed 

beams radiated by pulsed charge-current distributions on T>, still spinning 
at the uniform rate c/a. 
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1 The inertia density of an electromagnetic field 

The energy density u and the Poynting vector S of an electromagnetic field 
(E, B) in vacuum are given by 

u(x) = ^{E{xf +B(xf) and S(x) = E{x) X B(x), (1) 



where x = (x,t) € R 4 are the spacetime variables. We are using Heaviside- 
Lorentz units (eo = (mq = 1). S satisfies the vector identity 

S 2 = E 2 B 2 - (E-B) 2 , 

which implies that 

u 2 -S 2 = ^(E 2 -B 2 ) 2 + (E-B) 2 >0. (2) 

The electromagnetic momentum density is given by 1 

g(x) = S(x)/c, 

where c is the speed of light. For a relativistic particle with energy E and 
momentum p, the mass and velocity are given by 

2 

TO 2 c 4 = E 2 - c 2 p 2 > and v=°^ ^> v 2 < c 2 . (3) 

E 

1 In SI units [J99, page 261], g = S/c 2 . 
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Thus it makes sense to define the electromagnetic inertia density 



Z(x) = cr 2 ^u{xf - Jgixf = ^(E 2 - B 2 Y + A(E ■ B) 



and the electromagnetic energy Row velocity 



cS{x) 2 ^ 2 

v(x) = — tV, v(x) < c . 
u(x) 



(4) 



(5) 



T(x) and v(x) are local spacetime fields. Since E 2 — B 2 and E ■ B are Lorentz- 
invariant, 1{x) is a scalar field. 

From here on we set 



c=l, 

although c will be inserted into expressions occasionally to clarify their physical 
contents. 

The energy, momentum and inertia densities of a field can be expressed suc- 
cinctly in terms of the complex vector fields 



as follows: 



2u = F* ■ F = \F\ 



F(x) = E(x)+iB(x) 



2iS = F* x F 



21 = \F A 



(6) 



(7) 



The combinations E ± iB have been called Riemann-Silberstein vectors [B96] 
and Faraday vectors [B99]. They have been rediscovered many times and were 
explored extensively by Bateman [B15]. See also [K94, Chapter 9] and [K4]. 

By (4), 



2(x) = E(x) 2 - B(x) 2 = E(x) ■ B(x) = <^> F 2 = 0. 



(8) 



An electromagnetic field with F(x) — is said to be null at x. Nullity is a 
local, Lorentz- invariant concept. It is the Held counterpart of a massless particle. 
Indeed, 



\v(x)\ = 1 O X(x) = 0. 



(9) 



Electromagnetic energy propagates at the speed of light only at events x where 
the field is null. Elsewhere, it travels at speeds less than c and has a positive 
inertia density. 

Although this simple fact should be widely known in classical electrodynamics, 
I've been unable to find any clear reference to it in the mainstream literature 
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and in discussions with several knowledgeable colleagues. The sole exception, 
to my knowledge, is a brief note in [B15, page 6]. 2 

How are we to understand that electromagnetic waves, which are the ingredients 
of light itself, do not generally travel at the speed of light? The explanation is 
simple. 

• It is well known that every electromagnetic field becomes asymptotically null 
in the far zone, where its wavefronts become asymptotic to plane waves. But 
fields are generally not null in the near zone, close to their sources. Hence a 
generic electromagnetic field (i.e., one not chosen in a special way) has a non- 
vanishing inertia density in the near zone which decays to zero in the far zone. 
This explains why observed light (which is usually seen in the far zone due to 
its very short wavelength) is seen to propagate at a speed very close to c. 

• The inertia in the near zone can be understood as the result of interference 
between different parts of the field. As illustrated in Figure 1, parts of the near 
field propagate away from the source while other parts propagate back toward 
the source. That causes reactive energy to circulate back and forth between the 
source and the near field, giving rise to inertia. This has important consequences 
in antenna theory, where electromagnetic fields tend to loiter near the antenna 
to a larger or lesser extent. Antennas generating a great deal of reactive energy 
are inefficient radiators since the associated inertia tends to slow down their 
transmission time, effectively limiting the bandwidth. For a careful analysis of 
the reactive energy in the near zone of antennas, see [Y96]. See also [Wik-NF] 
for a less rigorous but easily comprehensible description. 

• An important aspect of nullity was discovered by Iwo Bialynicki-Birula [B3] . 
A generic electromagnetic field in free space is null along a set of 2-dimensional 
hypersurfaces S in spacetime since F(x, t) 2 = imposes two real conditions 
on the four spacetime variables (x,t). The time-slices St of S are (generically) 
curves in space, and these curves evolve with t. Bialynicki-Birula showed that 
the non-null field surrounding such a curve circulates around it, forming a vortex 
line along St- The above argument shows that when a field is not null in an 
extended region of spacetime, it travels at the speed of light only along such 
electromagnetic vortices. Elsewhere it travels at speeds less than c, although 
the energy flow speed everywhere approaches c in the far zone. 

• The characteristic property of fields that are null in an extended region of 
spacetime (not merely on isolated vortex lines, as in the case of a generic field) 
is that they are coherent, i.e., their various parts travel in unison and without 
interference. They are rewarded by having zero inertia and propagating at the 
exact speed c within the region of coherence. 

Various classes of globally null electromagnetic fields (where F 2 vanishes almost 

2 In quantum electrodynamics it is well known that exchanged photons are virtual (off the 
zero-mass shell in Fourier space), hence they can have any speed - possibly greater than c. It 
should be interesting to connect these classical and quantum ideas. They are in some sense 
dual: one is local in 4D spacetime, and the other is local in 4D Fourier space. 
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Figure 1: In the near zone, the energy in a generic wave propagates back to the source 
as well as away from it. This incoherence leads to reactive energy and an energy flow 
speed v < c, hence inertia. In the far zone, only the outgoing waves survive and the 
field becomes asymptotically coherent. 



everywhere, with the possible exception of singularities whose supports have 
zero measure) are known; for examples, see [B15, R61, T62, RT64]. However, 
no solutions with extended sources appear to be known which radiate null fields 
everywhere outside the source region. 

In this paper we construct a class of localized pulsed-beam fields whose charge- 
current densities live on a disk in space and which are null everywhere outside 
the disk. 3 We call these fields coherent electromagnetic wavelets. If the singular 
sources radiating such wavelets can be approximated in practice, they should 
be very efficient radiators. 



2 The scalar pulsed-beam wavelets ^(x,t) 

We now summarize some results regarding pulsed-beam wavelet solutions of the 
scalar wave equation. Detailed explanations can be found in [K3, HK9]. 

Consider the wave ^(x, t) radiated by the point source 4tt6(x) driven by a pulse 
go(t), which satisfies the wave equation 

nV(x,t) =4TTg (t)5(x), whcreD = 9 t 2 -A (10) 

is the wave operator. The unique causal (retarded) solution is 

*(x,t) = 9o(t - r) ^ 9o{t - r/c) . (II) 
r r 

The corresponding scalar pulsed-beam wavelet ^> is an extension of 'J to complex 
spacetime, obtained formally by displacing the point source from the origin to 

3 More accurately, the fields are also singular on the line forming the axis of the disk, which 
acts as a vortex jet. See section 9. 
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an imaginary point ia. It is defined by 

*(x,t) = ^^ (12) 

where t — t — is is a complex time, 

£ = y/(x — ia) 2 = V r 2 — a 2 — 2ia ■ x, r = \x\, a = \a\ (13) 

is the complex distance 4 from the imaginary source point ia to the real obser- 
vation point x, and g is the analytic signal or positive-frequency part of the 
driving signal go(t), i.e., 

1 poo i poo 

g(T) = — foe-™ fld(w)=— / due-^e- u 'g (u). (14) 

Here go (u>) is the usual Fourier transform. If <?o is square-integrable, then g 
is analytic in the lower-half complex time plane s > due to the smoothing 
factor e _ " s in (14). If go(u) decays rapidly as ui — > oo, then g(t — is) may also 
extend analytically into the upper half-plane. For example, in [HK9] we used 
the Gaussian pulse of duration d > 0, 

so that 

d^O => -> til. 

r 

This is the retarded wave propagator, a fundamental solution of the wave equa- 
tion: 

□#(x,t) = 47T(5(t)(5(a;). 

The associated pulsed beams ^(x,t) were shown to form a 'basis' for all radi- 
ated solutions of the wave equation by extending Huygens' principle to complex 
spacetime. Thus "J are 'wavelets' in two different senses: 

4 The idea of transforming a spherically symmetric solution of a partial differential equation 
with a point source to a cylindrically symmetric one by an imaginary displacement of the 
point source, and the ensuing complex distance, was introduced by Ted Newman in general 
relativity in 1965. Newman and his collaborators used it to derive the Kerr solution (neutral 
spinning black hole) and the Kerr-Newman solution (charged spinning black hole) from the 
spherically symmetric Schwarzschild and ReissnerNordstrom solutions, respectively [NJ65, 
N65]. Deschamps discovered the idea independently in 1971 in the context of the Helmholtz 
equation, where he use it to derive the time-harmonic complex-source beams [D71] from the 
spherically symmetric Green function, as in (25): 

> . 

c 

Probably the earliest use of £ was in 1897, when Appell used it to study the analytically 
continued Newtonian potential 1/f as a harmonic function [A87]. 
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• From a classical point of view, VP and its spacetime translates ^>(x — xq) are 
pulsed-beam versions of Huygens wavelets. The analytic continuation to 
complex spacetime breaks the spherical symmetry of the latter, deforming 
them to pulsed beams. 

• From a modern point of view, they are a family of localized functions 
paramaterized by their centers and widths which can be used to expand 
a space of functions, as in time-scale analysis [K94]. In this case, the 
centers are the real spacetime translation parameters xq, the widths are 
the imaginary translation parameters (a, s), and the space of functions 
represented as superpositions of \t are radiated waves. 

Returning to the general wavelet * in (12), we define its source q by 

&-KQ{x,t) = n$(x,t), (16) 

where the wave operator □ must act in the distributional sense to accommodate 
the singularities of VP. To understand g, note that the complex square root £ 
has branch points whenever its argument belongs to the negative real axis: 

r 2 — a 2 — 2ia ■ x < <^ r < a and a ■ x = 0. 

Since we are keeping the imaginary spacetime variables (a, s) fixed, the branch 
points of ( form the disk in M 3 given by 

V = {x : r < a, a ■ x = 0} = {x : Re C = 0}. (17) 

The formal displacement of the point source from the origin to ia transforms 
the point singularity of to a disk singularity of ^> on T>. \P is infinite on the 
boundary 

C = dV = {x : r = a, a ■ x = 0} = {x : ( = 0} (18) 

and discontinuous across the interior of T>, where Im £ has a jump discontinuity. 
We choose the branch of £ with Re C > 0, so that £ is an extension of the positive 
Euclidean distance, i.e., C — > T as a — s- 0. It then follows that \P is an extension 
of \& in the sense that (a) its source £>(x,t) vanishes whenever x £ T>, and (b) 

a^O, s->-0 => 2 Re V(x, t) -> - ^ = V(x, t). 

r 

Taking the real part is necessary to restore the negative-frequency component 
of go. The distribution q in (16) has been computed rigorously in [K0, K4a]. 

As its name indicates, ^(a;,^ is a pulsed beam radiated by the disk T> and 
propagating in the direction of a. To see this, write 

C = e-i»7, £>0. (19) 
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Choosing a cylindrical coordinate system (p, <fr, z) with z-axis parallel to a, we 
have 



hence 



and 



C 2 = (x - iaf = r 2 - a 2 - 2ia ■ x = f - rj 2 - 2i£r], 



a 2 = £ 2 — rj 2 and az = £77 



2 2 22 22 2f 2 , t-2 2\ <-2 2 

a p — a r — a z — a (a + £ — 77 J — 5 77 ■ 
Therefore (£, 77) are related to (p, z) as follows: 

a 2 p 2 = (a 2 + £ 2 )(a 2 - rj 2 ) and az = £77 



(20) 



In particular, it follows that 



• < £ < r and —a < rj < a, 

• the level surfaces of £ are the oblate spheroids 



P 



z 



= 1 



0?+^ £ 2 

• the level surfaces of rj are the semi-hyperboloids 



£>0, 



(21) 



H n = <, x 



P _ z_ 
a 2 — rj 2 rj 2 



= 1, Sgnz = Sgnr? L < t/ 2 < a 2 . (22) 



The far zone consists of points x with r»a, where the observer is far from the 
disk. By (13) and (19), 



a => £wr and rj ps a • x = a cos 0, x = x/r 



(23) 



Hence £ is a deformation of r and 77 is a deformation of acos#. For £ ^> a, 
becomes asymptotic to the sphere of radius r = £ and "H,, becomes asymptotic 
to the cone cos 6 — rj/a. The circle C (18) is the common focal set of all the 
spheroids and hyperboloids. See Figure 2. 

To see that \P(x, t) is a pulsed beam propagating in the direction of a, let go(t) 
be a pulse peaking around t — 0. Then at a fixed observation point a;, 



fl(r-C) 



1 

2^ 



(24) 



is a pulse peaking around t = £(x). Hence the wavefrants of ^(a^i) consist of 
the oblate spheroids 0%, just as the wavefronts of *f>(x, t) in (11) are the spheres 
r = t. 
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Figure 2: The real and imaginary parts of £ = £ — 177, together with the azimuthal 
angle <j>, form an oblate spheroidal coordinate system (£,77, 0) in R 3 . The above plot 
shows cut-away views of an oblate spheroid with £ = 0.7a, a semi-hyperboloid 7-L v 
with 77 = 0.8a (2 > 0) and another with 77 = —0.5a (z < 0). Also shown is the focal 
circle C with radius a at z = 0, whose interior is the branch disk T> of £. 



But the pulse amplitude is not uniform on 0^. Due to the factor e^^^^ in 
(24), it is strongest along the positive z-axis, where 77 = a, and weakest along 
the negative z-axis, where 77 = —a. In fact, if we follow a propagating wavefront 
along any given semi-hyperboloid H v (i.e., fix 77 = 770 and let £ increase as 
=* + &), then 

*(x,t) = — — — : , i»£o-«7o- 

The numerator remains constant, and the attenuation of \& is due entirely to the 
increasing distance, i.e., size of the denominator. This shows that \P propagates 
along H v , so each V. v is a diffraction surface for the pulsed beam. In the far 
zone, 0£ becomes a spherical wavefront (£ ~ r) and becomes a diffraction 
cone (77 sa a cos 0) for the beam. 

In the frequency domain, ^ is represented by 

% u (x) =&(w)^-=&(w)^-e*"> w >°> ( 25 ) 

where fe = w/c = w is the wave number and we have suppressed the factor e~ w;i , 
which controls the growth as lo —¥ 00. These are the original complex-source 
beams derived by Deschamps [D71]. In the far zone, 



r>a * w (x)~5o(w) e 

r 



&a cos 6 
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showing that the radiation pattern of the beam is 

^(0)=9o^)e kacos6 . 

Applications of complex-source beams in electrical engineering were developed 
extensively by Felsen and later extended to the time domain by Heyman et al. 
using analytic signals. See [F76, HF89]; also [HK9] and references therein. 

The spacetime analytic-signal transform, a tool for extending solutions of wave 
equations (including the Klein-Gordon and Dirac equations in quantum me- 
chanics) to complex spacetime and their interpretation as directed wave packets 
(relativistic coherent states), was introduced and developed in [K77, K78, K87, 
K90, K94, K3]. 

Remark. The analytic signal g(r) can be interpreted as follows. Whereas the 
point source 8{x) can be driven by an arbitrary time signal go{t), the disk T> 
must be driven by a separate time signal go{p, t) for each p € T>. However, these 
signals cannot be independent since the points p are rigidly connected to T>. By 
the magic of complex analysis, the single analytic function g(t — is) represents 
a coherent driving signal for the entire disk. The imaginary time parameter 
s represents a reaction time measuring the degree of correlation. Increasing s 
makes g(t — is) smoother, hence more correlated. In (24), for example, we must 
require that 

s — Tj > Va3, hence s > a = a/c, 

in order for the integral to converge when g is a general square-integrable signal. 
This implies that the correlations between all points on T> are related causally. 

3 The twisted null congruence of \I/ 

We now show that the oblate spheroids and semi-hyperboloids H.^ can be un- 
derstood geometrically as the wavefronts and diffraction surfaces, respectively, 
of 'light rays' (points moving in straight lines at speed c) emitted from the 
disk T> if the latter is assumed to spin at the extreme relativistic angular ve- 
locity c/a. This will later inspire the construction of coherent electromagnetic 
wavelets, which fully justify the interpretation in terms of light rays. 

We have shown that r H ri should be interpreted as a diffraction surface for the 
pulsed beam *f>, i.e., a surface along which = g(r — () is constant when the 
wavefront expands as £ = t. We now refine this interpretation. 

In Figure 2, imagine that H v consists of individual 'rays' propagating away from 
the source disk in the direction orthogonal to the wavefront . The location of 
the points tracing out the rays is given in oblate spheroidal coordinates (£, rj, (j)) 
with 

f = t, n = rj , <j>= 0o, 
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with 770 and 4>o fixed. We shall use £ as the time parameter to describe the 
motion. The trajectory of the ray in space is then 



and (20) gives 



Its velocity is 



where 



It follows that 



x o(Q = (p(O cos 0o,p(<Osin</> o ,z(£)), ( 26 ) 
a 2 P (0 2 = (e + « 2 )(« 2 - Vo), az(0 = Zr l0 . 
v (0 =d t x (0 = (p ? cos </>, p c sin 0, z ( ) 

Pt : = d s p= , z i = d i z = — . 

a l p a 



-1 = 4+4 = $^- (27) 

Since p(0) 2 = a 2 — T}q, the rays starting from the center of the disk (p(0) = 0) 
have r]o = ±a and travel in straight lines with speed |t?o| = 1- But all other rays, 
starting with < p(0) < a, have t]q = a 2 — p(0) < a 2 , hence their propagation 
speed is \vo\ < 1 and their paths are hyperbolic. 

Ideally, we would like our 'rays' to travel like light in free space: in straight 
lines, and at speed c. This can be arranged by allowing <fr to vary in time. Let 
the ray have the trajectory 

x(0 = (p(Ocos0(O,p(0 sintf(0,*(0), 
so that its velocity is 

v(£) = d t x(0 = v (Z) + p(04>'(04> 

where 

4> = d^p = dff, (cos</>, sin (j>,0) = (— sin</>, cos(£, 0). 
We require that \v(£)\ = 1, hence (27) gives 



P 2 r = 1 - 



e+a 2 e + a 2 r {e+a 2 f 

That is, the points of the circle 
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propagate precisely at the speed c = 1 if and only if spins at the angular 

velocity 



w± (O = 0'(O=±tfrr^ = ± 



£ 2 + a 2 £ 2 + a 2 



(28) 



Furthermore, since </>'(£) is independent of 770, a/i the circles Q-i ))0 with given 
£ spin at the same angular velocity and the oblate spheroid spins as a rigid 
body. We say that the spheroid with angular velocity lo + has positive helicity 
and that with angular velocity w_ has negative helicity. 

Since the branch disk T> is the limit of as £ —> 0, it follows that V spins 
rigidly at w±(0) = ±c/a. This is precisely the angular velocity at which the 
points on the boundary dT> move at the speed of light! More will be said about 
this in Section 4. 

Let us now compute the trajectory of a single point of Q_.^ . We can set 4>q = 
without loss of generality. Integrating (28) gives 

0(O = ±tan- 1 (C/a), (29) 

hence 

cos 0(0 = ° and sin 0(0 = ± 

V£ 2 + a 2 V£ +« 



Furthermore, ap(£) = y' (£ 2 + a 2 ) (a 2 — ?y 2 ), so 

p(Ocos0(O = ,Ja?-T$=p , p(Osin0(O = ±^, *(0 = 



The trajectory of the ray starting at the point on T> with p = p and = is 
therefore 

x(0 = L, ± — , — ) - x(Q) + «±£, (30) 
\ a a J 

where the velocity is given in Cartesian coordinates by 

v± = a -1 (0,±p Q ,i7o) , v 2 - = 1. 

Since the choice 0(0) = is arbitrary, we have proved the following result. 



Theorem 1 Every point emanating from T> which follows the exploding wave- 
front 0£ with £ = ct as it spins at the angular velocity w±(£) while constrained 
to Una moves in a straight line at speed c. 

The trajectories of such points might be interpreted as 'light rays,' except that 
they are associated with a solution of the scalar wave equation rather than 
Maxwell's equations. This will be corrected later. 
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The existence of two families of straight lines generating H v is a well-known fact 
of geometry: every circular hyperboloid of one sheet is a doubly ruled surface. 
Our semi-hyperboloids are ruled by the rays, which are half-lines emanating 
from a circle on T> as shown in Figure 3. The rays move upwards if 77 > and 
downwards if 77 < 0. 5 The motion is vertical if 77 = ±a (p = 0). It acquires 
a horizontal component when < \r]\ < a, and becomes purely horizontal as 
37 — ^ (p a). For r\ = 0, the trajectory begins on the branch circle C and is 
tangent to C, remaining in the xy-plane. That is so because C is already moving 
at the speed of light, so the vertical velocity component must vanish. The union 
of all such horizontal rays forms the degenerate semi-hyperboloid 

Ho = {x : Tj = 0} = {x : r > a,z = 0}. (31) 





Figure 3: The semi-hyperboloids H v are doubly ruled surfaces, with the straight lines 
interpreted as 'light rays' radiated by T> and traveling at speed c. 



Let us now compute the velocity vector of the ray going through an arbitrary 
point x £ K 3 . Identifying t again with £, define 

dx dx dp dx dz dx d<fi 



u(x) = 



Using (20) and 

dx 
dp 

we obtain 



dx „ dx 

dz J dd> 



p4>, 




(32) 



In terms of the orthonormal oblate spheroidal basis (£, 17, 0), this is seen to be 

- ! e 



u±{x) = £4 



±01 



(33) 



5 In our case, the ruling for r) < has the opposite helicity from that for r] > 0. That is 
because when the disk spins counterclockwise as seen from the top, it spins clockwise as seen 
from the bottom. 
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As expected, the coefficient of fj vanishes since the rays propagate with constant 
77. Furthermore, u± is not orthogonal to the wavefront since 



£ a 2 — n 2 

Note that the ^-component of u± does not depend on the sign of 77. This con- 
firms that the upper and lower semi-hyperboloids (±77 > 0) are ruled oppositely. 
Both (32) and (33) show that \u± \ = 1, as expected. 

The above picture of radiation from a rotating disk explains both the wavefronts 
0(: and the diffraction surfaces H v . With some poetic license, let us call the 
emitted particles 'photons' even though "J is a scalar potential and not an elec- 
tromagnetic potential. 6 As shown above, the envelope of all the rays emanating 
from a fixed circle of radius po < a in T> is the pair of semi-hyperboloids / H± riQ 
with 770 = \J a 2 — Pq . However, if we take a snapshot of all the photons at fixed 
time t — £ > 0, we find that they form the wavefront O^. This is illustrated in 
Figure 4. The radiation from a spinning disk thus explains the level surfaces of 
both £ and 77, showing the complex distance C = £ — if) is a natural tool for its 
analysis. 




Figure 4: This picture gives a geometric explanation of H v and in terms of the disk 
T> spinning at the extreme relativistic angular velocity c/a. Left: The envelope of the 
rays emitted from a circle or radius < p < a on T) is the diffraction semi-hyperboloid 
H v with 77 = ±\/a 2 — p 2 . (We show only T-L v with 77 = 0.7.) The larger p, the greater 
the horizontal speed and the more twisted is the bundle of rays generating H v . Right: 
The 'photons' emitted from the entire disk T> at t — form the wavefront at t — £. 



It can be shown that the vorticity of u± is given by 



277 

"£ 2 + 77 2 



u± = V x u± = ± - P) ' 2 u ± . (34) 



6 In Section 5 we shall reinterpret ^7 as the scalar component of an electromagnetic 4- vector 
potential, so that the terminology, while still classical, is not entirely inappropriate. 
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Thus u± is a Beltrami Geld, curling around its own axis with hclicity ±. 

There is a natural four-dimensional (spacetime) formulation of the above. Con- 
sider the two stationary 4- vector fields u±(x) with contravariant and covariant 
components 

u£(aO = (ti±(x),l) u±(x) = (-u ± (x),l), M= 1,2,3,4. (35) 

They are nul2, i.e., their Lorentzian squares with respect to the Minkowski 
metric rj^ = diag (—1, — 1, — 1, 1) vanish: 

2 v U ± M -i 2 r\ 

U ± = TJ^Uj-Ug = U± = 1 - «± = 0, 

where we sum as usual over the repeated index /i. u± are the 4- velocities of the 
twisted ray bundles with helicity ±, and they determine the rays by 7 

x^{x,t) = (x,0) +tu^(x) = (x + tu ± {x),t). (36) 

In general relativity, ray bundles (possibly in curved spacetime) such as (36) are 
called null congruences, and they play an important role in the construction and 
analysis of nontrivial solutions to the Einstein and Einstein-Maxwell equations; 
see Section 11. Note that our rays are not orthogonal to the wavefronts O^. This 
is a characteristic of twisted congruences, which are characterized byM-Vxtt^ 
0. 8 In our case, 

u ± • V x u ± = ± ^ V ^ Vz ^ 0. (37) 

4 + Tf 

The congruence is twisting everywhere outside the xy-plane, where its rays span 
the degenerate hyperboloid Ho (31) of half-lines tangent to C = &D. 

Remark. It will be shown in Section 11 that the null congruence in (35) is 
identical to the null congruences of the Kerr and Kerr-Newman metrics in gen- 
eral relativity [R61, T62, RT64, NJ65, N65]. This is an intriguing fact whose 
significance remains to be fully understood since the Kerr metric and electro- 
magnetic fields are stationary while VP is time-dependent. 



4 Newman's 'magic' electromagnetic field 

In 1973, Ted Newman published a paper [N73] detailing the electromagnetic 
part of the Kerr-Newman solution of the Einstein-Maxwell equations in general 
relativity [N65], which represents a spinning, charged black hole. Although the 
Kerr-Newman solution is quite complicated, its purely electromagnetic part is 

7 Here t ranges over all real numbers and there is a unique to for which x + t$u± £ T>. For 
t > to, x±(x,t) has helicity ±, while for t < to it has helicity =p. 

8 In general, a four-dimensional version of (37) must hold. But since u± depends only on 
x and the time component u l ^_ = 1, this reduces to (37). 
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a simple static field in flat spacetime. Newman begins with the analytically 

continued Coulomb field obtained by replacing the Coulomb potential r _1 by 
-l. 



c 



E(x) = -V± = (38) 



He then notes that E has the following multipole structure: 

Re E = (electric monopole) + (electric quadrupole) + • • ■ (39) 
Im E = (magnetic dipole) + (magnetic octupole) + • • • . 

For example, using the far-zone approximation (23) gives 

r > a => E = ^ - — + i + 0(a 2 /r 4 ). 

(r — ix • a) 13 r J r 6 

The first term is a real electric monopole of charge 5 = 1, and the second term 
is an imaginary magnetic dipole of magnetic moment a. Newman therefore 
interprets E by writing it as 

E(x) = E(x) +iB(x), (40) 

where E a real electrostatic field and B is a real magnetostatic field. 

In [K4b], I computed the charge- and current densities of the field (E,B) as 
the distributions in M 3 defined by 

Ang(x) = V • E(x), 4ttJ{x) = -iV xE(x) - dt^x) = -iV x~E(x), (41) 

which is just the complex form of the usual Maxwell equations. The use of 
distribution theory is necessary because E is discontinuous on the branch cut 
T> of C and singular on the focal circle C — dT>, where £ = 0. I showed that: 

• g and J are real distributions supported on I?, so that the homogeneous 
Maxwell equations are satisfied and thus no magnetic charges or currents are 
introduced by the complexification. 

• The surface charge-current density on the interior of T> is given by 

2 ° « = ~rj4^- (42) 



( a 2_ /0 2)3/2' ( a 2_ /9 2)3/2- 

This is the charge-current density of a rigidly spinning disk rotating in the 
counterclockwise direction cj) with angular velocity 9 

«=W = 1 = C . (43) 
up a a 



9 The reader may wonder how c has entered an expression obtained by analytically con- 
tinuing the Coulomb field. The answer is that the dimensionally correct expression for the 
current (41) is 

AnJ(x) = -icV x E(as) - d t E(x) = -icV x E(sc), 
so K must be multiplied by c. 
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The focal circle C therefore spins at the speed of light. 10 Note that the speed of 
a point on T> is v = cp/a, hence (42) can be rewritten in a form displaying its 
relativistic nature: 

K = ^ (44) 



a 2 (l-v 2 /c 2 ) 3 / 2 ' a 3 (l -v 2 /c 2 ) 3 / 2 

The surface sources (42) can be obtained easily from the discontinuities of E 
and B across the branch cut T>. As x — > T> from above and below, we have 

p < a, z — > 0* =>■ C — ^ T»v ft2 — P 2 , x^pp=p. 

The upper and lower boundary values of E are thus 

E± _» — — 

± 2 (a 2 - i o 2 ) 3 / 2 ' 

showing that 

± _ t 01 D± 



( a 2_ p 2)3/2> -r ( a 2 _ p 2)3/2 ' 

and the discontinuities are 

8E = E + - E 



2a 



5B = B + B 



(a 2 - p 2 ) 3 / 2 
2p 



(a 2 - p 2 ) 3 / 2 ■ 

The jump conditions at an interface [J99, page 18] thus give 

proving (42). Furthermore, 

z • SB = and z x 5E = 0, (46) 

proving that the magnetic charge-current density on T> vanishes as required by 
the homogeneous Maxwell equations. 

However, these surface sources do not give the entire story. The jump conditions 
do not account for a 1-dimensional charge-current distribution along the bound- 
ary C = dT> since the usual method [J99, page 18] of infinitesimal pillboxes and 
rectangular line integrals fails there. This failure is reflected by the fact that 
the total charge of a is q = — oo, whereas the field E was derived under the 
assumption that q = 1. 

10 This bizarre property is related to the fact that E is the flat-spacctime limit of the electro- 
magnetic field around a spinning charged black hole. That also explains the singular nature 
of q and K on C. See [N65, N73]. 
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The solution to this puzzle is to compute the complete volume charge-current 
densities, defined as distributions using (41). When they are computed rigor- 
ously (this requires a regularization) , we recover (42) on the interior of T> but 
also find that &D carries a linear charge-current distribution of charge +00, 
precisely so that the total charge remains q = 1. 

It is interesting to compute the inertia density I and energy flow velocity v of 
the Newman field E. By (7), 

^ = ^' henCeI= 2^=2FTW' ^ 

which is infinite on the focal circle C where £ = 0. Its energy density and 
Poynting vector are 

2 1 1 2|C| 4 2(£ 2 + 77 2 ) 2 2(£ 2 +?? 2 ) 2 1 ; 

a x x an - 



so the energy flow velocity is 



S - 2ap 4>. (49) 



U 2a 2 + £ 2 — rj 2 
It follows from (20) with X = a 2 + £ 2 and Y = a 2 - rj 2 that 



v 2 = '- = 1 - — — < 1 = c 

(X + Y) 2 (X + Y) 2 ~ ' 

with equality only when £ = 77 = 0, i.e., x £ C. 

Conclusions: 

• Note the similarity between the charge flow (43) on V and the twisted con- 
gruence derived from ^ in Section 3. In both cases, T> is seen to spin at the 
extreme relativistic rate u> — c/a. This will be the inspiration for reinterpreting 
VP as the scalar part of an electromagnetic 4-potential (A, whose field is a 
pulsed version of the Newman field E (Section 5). 

• The energy field circulates about the z-axis without propagating out (49), 
which is expected since the field is non-radiating. Its angular velocity at x is 

n(x) = - = - „ 2 ° — (50) 

y J p 2a 2 +£ 2 -n 2 y ' 

For fixed £ > i.e., x € Of, il(x) is a minimum on the equator (rj = 0) and 
increases steadily to its maximum at the north and south poles as \r)\ increases 
from to a. Hence the energy flow forms a vortex along the z-axis. 



AXY (X - Y) 



2 
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• On T> (£ = 0), we have 

. 2a 2a 2ac 

"IP) = o = "o o = ~n o) P < a. 

2a 2 -if a 2 + p 2 a 2 + p 2 

The angular velocity of the energy flow on T> differs from the uniform angular 
velocity u) (43) of the charge flow. It increases from its minimum value f2(a) = 
c/a = on C to its maximum value 51(0) = 2lj at the origin, which is the center 
of the vortex. 

• The real Coulomb field is stationary, as expected: 

a -> =*> v ->• 0. 

The rotation thus comes directly from the analytic continuation r — > £. 



5 Prom scalar to electromagnetic wavelets 

In this section we introduce the idea of electromagnetic pulsed-beam wavelets 
as time-dependent (pulsed) versions of Newman's static field E(x). Like the 
latter, they will be formulated in complex form as 

F(x,t)=E(x,t)+iB(x,t). (51) 

We begin by interpreting the pulsed beam ^(x,t) of Section 2 as a scalar poten- 
tial and computing a complementary vector potential A(x,t). Both and A 
are complex. The real field is then obtained as follows. First define the complex 
fields by complexifying the usual expressions, 

E = - W -d t A B = V x A, (52) 

and let 

F = E + iB. (53) 

Then (51) implies 

E = Re E - Im B B = Re B + Im E. 

Since ^> in (12) is retarded, we require that the four-vector (A, ^) satisfy the 
Lorenz gauge condition 11 

V-A + d f * = 0. (54) 

11 Other gauge conditions do not respect causality. While this is acceptable for potentials 
due to the gauge invariance of the fields, it is not desirable in our case. For example, the fact 
that ^ depends on the complex retarded time r — £ gave us both the wavefronts and the 
diffraction semi-hyperboloids H v - 
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This relates the unknown vector potential to the known scalar potential. We 
furthermore require that the current density must vanish outside the branch 
disk V of (, which is also the support of the charge density g defined in (16): 



AttJ{x, t) = DA(x, t) = Mx^V. 



(55) 



Since we expect A to have the same time dependence as given by the factor 
g(r — C), we make the Ansatz that 



A(x,t) = ^(x,t)w(x) 



(56) 



for some static vector field w(x). This will now be used to solve Equations (54) 
and (55). Since 



w = - g ' (r ~ c) v C - 



£2 s 



9_ , 9_ 

c c 2 



VC, 



(57) 



we have 



V • A = V* • w + W • w = - ( ~ + ) VC • w + ^V • w. 



Hence (54) requires 



9_ 9_ 

c c 2 



VC • w + ^V • w + j = 0. 



Because this must hold for all analytic signals g, we have 



V • A + d t ^ = 0Vg <f> VC • w = 1 and V • w = 



(58) 



Using (56) in (55) gives 

□ = (□§)«; - 2(V* ■ W)w - i>Aw = 0Vx(£V. 
Since = outside T>, the first term vanishes and (57) gives 



£+9 

V C C 2 

Again, this must hold for all g, hence 



2 ( ^ + ^ ) (VC • V)w - ^Aw = Vx i V. 



fJ 

r 



(VC • V)w = and Ato = Mx (£ V. 



(59) 



Finding A thus reduces to finding a vector field w(x) satisfying (58) and (59). 
Amazingly, although this system appears to be overdetermined, we shall find a 
reasonable set of solutions w. 

Note: If w = w(x, t), then (58) remains the same but (59) generalizes to 



d t w + (VC • V)«? = and Dw + -d t w = 0. 



(60) 
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6 The complex spheroidal frame defined by ( 



Our main tool for computing w(x) will be a complex moving frame [F63] in C 
which conforms entirely to the spheroidal geometry defined by £. Let 

- x — ia pp + zz — iaz, pp + zz 

C = VC = = = , where z = z — ia. 

c c c 

This is the complexification of the real unit vector f = Vr, and it satisfies 



Note that 



{Vi-iV-qf = 1 =*> |V^| 2 - |V?y| 2 = 1 and V£-Vn = 0, 

confirming that the levelsurfaces and H v are orthogonal. Define the com- 
plexification ■& of the polar angle 9 by 



sin#=~, cosi? = -, sin 2 d + cos 2 d = 1, 

c c 



so that 



C = P sin $ + z cos 



Next, define the complexification of the unit vector 9 = cf) x £ by 



i) = cf> x C = p cos # — z sin i? 



(61) 



(62) 



(63) 



whose real and imaginary parts can be shown to be given by 

- = + anp j]r§ - ajp 

e+v 2 e+v 2 ' 

The three vectors 

Ui = C, u 2 = 1?, u 3 = 4> 
satisfy the bilinear orthonormality conditions 12 

u k ■ ui = Ski 

and form a right-handed complex-orthonormal frame in 
w(x) has the unique expansion 



w = ^2{w ■ u k )u k . 



(64) 
(65) 

Every vector field 
(66) 



12 Not the sesquilincar orthonormality condition ut ■ ui = S k i based on a hermitian inner 
product like the one used in quantum mechanics. 
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The entire basis is determined by the complex distance (, and thus it conforms 
completely to the associated complex spheroidal geometry. The use of in 
defining # does not spoil this because ( is independent of 4>. Consequently, 
the electromagnetic field (51), including its polarization, will conform to this 
geometry. As we shall see, that makes it possible to construct coherent fields 
with vanishing inertia because the polarizations on nearby rays do not clash. 

Let the unknown vector field w have the expansion 

w(x) = a(x)C + ^(x)^ + j(x)4>. (67) 

We want to use Equations (58) and (59) to compute (a, f3, 7). The most difficult 
one of these would appear to be the first equation in (59), i.e., 

D ( w = where D c = C ■ V, (68) 

since the directional derivative generally acts on the basis vectors as well as 
their coefficients. Note that the gradient operator can be expressed by 

v = vc o c + vc*a c * + a = <d c + c*a c * + 1 (69) 

where 

= + dl = ~{d t -id n ) (70) 

are the complex derivatives. Hence (68) becomes 

D ( w = where D c = C ■ V = d c + |C| 2 9^. 



A glimmer of hope appears when we consider the real limit: 

a -)• => D r f • V = d r . 

Since the spherical basis vectors do not depend on r, we have 

d r r = d r 6 = d r 4> = 0, (72) 

and d r does indeed differentiate only the coefficients of a vector field expressed 
in that basis. 

The complex case is more subtle. While the real variables r and 9 are indepen- 
dent, C and •& cannot be independent as complex variables since they arc both 
functions of (C,??)- 13 But somewhat miraculously, it turns out that the essential 
property (72) survives complexification. 

13 The precise relation between f and i? is 

4aCsintf = 4ap = ^4a 2 +4£ 2 v / 4a 2 - 4r? 2 = Ma 2 + (C + C*) 2 J^a 2 + (C - C*) 2 - 



(71) 
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Theorem 2 (a) The complexified polar angle $ is constant with respect to dif- 
ferentiation along the complex direction £. That is, 



D c $ = 0. 



(73) 



(b) Consequently, the basis vectors {C,#, <A} are are also constant with respect 
to D c : 



D c £ = D c tf = D c cf> = 0. 



(74) 



Proof: (a) Differentiating £sin$ = p with respect to £* and C gives 

(d^-d cos ■& — d^p sin ■& + (d^tf cos A = d^p 

From (20) and (70) it follows that 



d* cP = 



< > r 2iap 2iap 
The first equation in (75) gives 



cos $ 



d c p = 



Hap 



% ziap 



and the second equation gives 
(Oqi) cos $ : 



(z* p C, 2 z* — Hap 



2iap C 2iapC, 



Inserting Q = p + z and —2ia = z* — z, this reduces to 



which proves that 



d c $ = |C| 2 ^77 = -|C| 2 <9^, hence D c ti = 0. 



To prove (b), note that 
by (62) and (63), hence 



2iap 



and 



D C < = d c C + K\ 2 o* c c = # + |C| 2 #t? 



0D c i9 = 0. 



y fS - u yuQU -r |S| 

Furthermore, since <f> is independent of £ and £*, we have 

D c = and D c d = D c (</> x C) = x = 0. 



(75) 



(76) 



(77) 
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7 Computing w(x) 



We are now ready to compute 

w(x) = a(x)C + f3(x)& + ~f(x)4> (78) 
by enforcing Equations (58) and (59), which are here summarized as 



(a) C- io = l, (b)V •«>=-, (c)D c w = 0, (d) Aw = 0. 



(79) 



• Equation (79) (a) immediately gives a(x) = 1. 

• By (71), D C C = 1, £> f C* = ICI 2 , and 

s 2%ap 2iap 

Of the four variables £, £*, d, only $ is constant under Hence Theorem 2 
shows that /3 and 7 are functions of the angles </>). That is, they are analytic 
in just as ^ is analytic in £: 

IB = C + 0(0, </>)# + 7(l?, 

• Assume that A, like \&, is axisymmetric, i.e., its components are independent 
of</>. Then 

w = C + + 7(^)0 (80) 

and 

D c w = (D c /3)3 + (D ( j)4> = P'{D Q $)$ + ry'(D c &)4> = 

as required by (79) (c). We now list some relations that will be needed in solving 
the remaining equations (79) (b) and (d). 



Proposition 1 



vc = C Vx( = o v-c = 7 AC -7 A C- , 



2 2 - 2C 

- AC=- AC = - — 



^ n # „ a <t> „ a cos A n cot A s 1? + sin 2d C 
Vtf = - V x = ^ V • 1? = A?9 = — — Ai9 : 



C C p C 2 p ; 

V(f>=- Vx<f>=- V-0 = O A0 = O A0 = — ^. 
P P P 
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Proofs of the less obvious results are as follows: 



£ cos $ Vi? — p — (p sin $ + z cos i?) sin d = p cos — C cos $ sin i? = i9 cos $ 



• AC = AVC = VAC = V> - , • 



px (A z 

• V x $ = V x (pV0) = p x = — - = - 

P P 

• V-i9 = V-(0xC)-C-Vx0= 

P 

• v • 4> = v • (pV^) = p-4> + pA(j> = o 

,2 _2C 

c c 2 

Thus 

2 - 2 /?cos?9 /?' 
V-io=-+ / 3V-tf + /?' W • + 7 V • + 7 ' W •</>=- + + 4r 

C C p C 

and Equation (79) (b) becomes 

/3cosi? _ 1 

To solve this, multiply by Csini? = p: 

/3cos$ + fi' sintf = d$((3su\$) = — sintf, 

hence 

/3sintf = cos + n (81) 

where K is constant of integration. Thus we have used Equations (79) (a)-(c) 
to find 

w = C + + 7 where (3 = cottf + kcsc$. (82) 
It remains only to enforce equation (79) (d). We have 
Aw = AC + (A/3)tf + 2(V/3 ■ V)# + /3Ai? + (A 7 )0 + 2(V 7 ■ V)0 + 7 A0 

= AC + (A/3)S + 2/3' (Vi? • V)# + /3Ai9 + (A 7 )0 + 2 7 '(Vtf • V)0 + 7 A0. 

We find 

ft = - esc 2 ■& - « cot d esc i?, A/3=^, 

P 2 



£ • I C3 ( + C a c + - d I = — d c # = 



But 



9^ = = 2^-& 
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hence 



Since (Vi? ■ V)</> = 0, we have 
2 

, C + 

P 



e * 1 p 2 c 2 p 2 



C 2 ' 



i? + Csin2i9 



A7 - 



0. 



By yet another miracle, the coefficients of £ and i) both vanish and we are left 
with 



Ato 



A 7 -^ 
P 2 



0. 



(83) 



Thus (79) (d) requires that 7 satisfy the same equation as j3, for which we have 
already found a general solution: 

Aw = =s> A7=^- =s> 7 = Acotz? + /xcsci?. (84) 

P_ 

We have proved the following result. 

Theorem 3 The general axisymmetric solution of Equations (79) is 

w = C + (cot-d + kcsc$)# + (Acottf + pcscd)4> (85) 
where n, A, /1 are arbitrary complex constants. 

Remark. The solution (85) is singular on the branch disk T> since £, hence also 
"d, is discontinuous there. In fact, 



=>■ cos $ -> ± 



a 



y'a 2 - p 2 



z — > —la 

sin d — > ±i 



P 

\Jo? - p 2 



In addition, w and thus A, is singular on the z-axis. To bring out this singu- 
larity, recall that sini? = hence (85) can be rewritten as 



- c - c 

w = C + - (cos + k)& + - (A cos # + n)<j>. 
P P 



(86) 



8 The general pulsed-beam fields F± 



According to (56) and (86), the general axisymmetric vector potential comple- 
menting \E" is 



A = -C+ -(cosi? + k)£+ ^(Acostf + /x)0, 
C P P 



(87) 
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where we have used the shorthand g = g(r — £) and k, A, /i are arbitrary complex 
constants. 

Let us compute the complex electric and magnetic fields defined in (53): 
E = - W - d t A 



9_ a_ 

c e 



C- ^rC- -(cos?9 + k)i9- ^(Acosi? + ^)<£- 

C p p 



Note that the radiating longitudinal component (<//C)C cancels, leaving only 
the non-radiating longitudinal component (<?/C 2 )C : 



E= ^C- -{cosi} + K )-d 
C P 



-(A COSZ? -(- jlt)0. 



(88) 



The radiation is seen to be confined to the (complex) transversal directions # 
and <j). The 1/p decay of the radiating terms indicates that the singularity along 
the z-axis carries a charge-current distribution and thus represents a filament 
attached to the branch disk. This will be investigated in detail elsewhere. 



For computing B = V x A, it will be convenient to rewrite A as 
A = + #(cot?9 + «csci?)Vi9 + g(Xcosd + /i)V</>, 



(89) 



where we have used i9 = (Vfl, 4> — and C/p — csct?. Then 

B = -g'(cottf + «csctf)C xVi? - ^'(Acosi? + /z)C x V0 - Xg sin^Vi? x V(/> 
= -^(cottf + Kcsctf)</>+ — (Acostf + u)tf- Aff Sm ^ C 
whose radiative structure is seen more clearly in the form 



B 



c 2 



C+ ~(Acostf + /i)i?- -(costf + 



P 



(90) 



Note the close similarities between E and B, down to the non-radiating longi- 
tudinal parts. These similarities combine to give 



F± = E ± iB = p± + (q± - P± cos 0) 9 -y± 

C p 



where 



p± = 1 =p iA, q± = — k ± z/i, (^-t = i9 ± icj). 



Note that 92+ and <J5_ are nuiJ vectors: 



<^± = ■&' - <p~ ± 2itf ■ (f> = 0. 



(91) 

(92) 
(93) 
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They form a complex basis of helicity ± for the the directions i9, <f> orthogonal to 
£, and will play an important role in the representation of our coherent wavelets. 
Note also that 



± iV0 = Vip± 



where 



hence 



p sin'i? 

V>± =lntan(i?/2)±i0 = ln (e ±# tan(i?/2)) , 



P 1( ^± ' da; 



sin 



± id<p = dip± . 



This can be used to write 



F± ■ da; = p± — d( + {q± - p ± cos d)g'di/>± . 



(94) 



(95) 



(96) 



(97) 



Remark 1. It might seem that using both complex fields F± in (91) is super- 
fluous. But since E and B are complex, F + and F_ are not related by complex 
conjugation, thus both fields are needed to recover E and B: 

2E = F++F_ 2iB = F + -F_. 

On the other hand, to compute reai fields, we need either F + or F_ but not 
both. Define the two real fields (E + ,B + ) and (E_,B_) by 

F± = E±±iB ±) 

so that 

E± = Re F ± = Re E =F Im B (98) 
B ± = ±ImF±= Re B ± Im E . 

Clearly these two fields are inequivalent, so F + and F give two independent 
real solutions to Maxwell's equations. Each of these solutions corresponds to 
a class of equivalent complex fields. Namely, the real field (E± 1 B±) remains 
unchanged when 

E^E + E' ± and B ->• B + B' ± where B' ± = ±iE', (99) 

since this leaves F± invariant: 

F± -> (E + E±) ± i(B ± iE'±) = E±iB = F±. 
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If the changes are induced by changes (A' ± , ^' ± ) of the potentials, then 

-V#± - d t A' ± = E' ± = =rB' ± = =rV x A±. (100) 

Taking the divergence of both sides and using the Lorenz condition gives 

-A#± - d t V • A' ± = -A*' ± + <9 t 2 #' ± = nStt'± = 0. (101) 

Similarly, taking the curl of (100) yields DA± = 0. The addition of a sourceless 
term (A' ± ,#' ± ) to the 4-potential (A, thus leaves the real field (E±,B±) 
invariant. However, this does change the field (E T , B T ) of opposite helicity: 

Fzp — > (E + E' ± ) =F «(B ± iE' ± ) =F T + 2E' ± 
=> E T -> E T +2 Re E', B T -> £ T T 2 Im E'. 

Furthermore, note that (.E±,.B±) remains unaffected outside its source region 
even if (Aj_, , i'j_) has a charge-current density whose support is contained in 
that of (A, ^). This is a new kind of 'gauge freedom' belonging to the complex 
representation (91). In Secton 11 we shall encounter an example of a 'pure 
gauge' field where A^O but F± = 0; see (120). 

Remark 2. Equation (91) shows that the radiating parts of the two indepen- 
dent fields (E+, B + ) and (E-, B-) associated with F + and F_ have helicity + 
and — , respectively. This confirms that both fields F± are needed, as explained 
above. 

9 Coherent electromagnetic wavelets 

The inertia density (4) of the field F± is 

z* - 

Since £ • <p± = and (p± — 0, (91) gives 

Fl=pi 9 -tL^f. (102) 

Thus one, but not both, of the fields F + and F_ can be made null: 

X=-i =*> p+=0 F + = q+ 9 ® Cp+ => F 2 + =0 (103) 

P 

X = i => p-=0 F = g _ - ^ - ^ y_ F^=0. 
Up to a complex scalar factor, we have thus found two unique null fields. 
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Theorem 4 The field F± is null if and only if X — =R> and then 



F±(x,t) = g/(T C) y ± , F ± • dx = </(r - C)dV± 



(104) 



or a complex scalar multiple thereof, where if>± is given by (95). 



Definition 1 We call F± the coherent electromagnetic wavelets with helicity 
± and driving signal g. 

Recall that A contained three independent complex parameters k, A, fi. To 
make F± null, we must choose A = =fi, respectively. That leaves two complex 
parameters, given now by q± in (103). To reconstruct the complex fields E and 
B, we need both the null field F± and the non-null field F T , hence both q + and 
q_. However, if we are interested only in null fields, then we need only q + (for 
positive helicity) or q_ (for negative helicity). 

Theorem 5 F± is an an electromagnetic pulsed beam propagating in the direc- 
tion of a, with wavefronts exploding as £ = ct, and diffracting along "H^. 

Proof: Apply the argument used with Equation (24) to interpret f as a scalar 
pulsed beam to the pulse-shape factor in (104): 



g'(r -Q = -± I dwwe-^HJe-^)^^. 

Z7T 



DO 







Remark. Since (p± = $ ± icf> and 

zp — pi 



-d = 



c 



by (63) and (61), F± is singular on the branch disk T> of £. Furthermore, (104) 
shows that it is singular along the z-axis p = 0. F± is analytic elsewhere, and 
this implies that its charge-current distribution is supported on T>, representing 
a spinning charged disk, as well as the z-axis, representing vortex jets generated 
by the spin of V in the positive and negative z-directions. However, as shown 
in Theorem 5, the jet along the negative z-axis is much weaker due to the 
pulse-shape factor g'(r — Q. 

The sources of F± will be studied elsewhere. 



10 The real and complex null congruences of F± 

The null fields (103) define both real and complex twisted null congruences in 
spacetime. The real congruence follow from the nullity condition 

(E ± ±tB ± ) 2 = S 2 ± = \E ± xB ± \ 2 =u\ = \(E\ +B 2 ± ) 2 , 
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which shows that the energy propagates at the speed of light: 

S± 



v± 



u± 



v\ = % = 1 = c 2 . 

Ill 



Theorem 6 The energy flow velocity for the coherent wavelets (103) is 



v±(x) = - 



£ a 2 — rj 2 „ 77 ^ a 2 — 1] 



a V £ 2 + a 2 



p+ -z± 



£ 2 + a 2 







(105) 



e + a 2 



e + a 2 



which is identical to the velocity field (32)— (33) associated with the scalar wavelet 
"J. The 4-vector field v± = (v±, 1) thus forms a twisting congruence in Minkowski 
spacetime. This can be expressed as a differential form: 



, ± (x).dx = |J±4de ± ^ 

£ J + gt a 



(106) 



z = z — la. 



Proof: Recall from (62) and (63) that 

CC = pp + zi £i9 = zp — pi 

It follows that 

|C| 2 i9* • S = r 2 + a 2 = |C| 2 |C? = e - V 2 + 2a 2 
|C| 2 i9* X 3 = 2iap4>. 
To simplify the notation, let 

X = £ 2 + a 2 and Y = a 2 - rj 2 

so that 

XY = a 2 p 2 X + Y = £ 2 -rj 2 + 2a 2 X - Y = £ 2 + n 2 = |C| 2 . 
From (p± = d ± icj> it follows that 
|C| 2 ^-^±=2X 

|C| 2( P± x <£± = i2t £,pp + (Xr]/a)z ± apcf> ■ 
For the null field F± = (g' / ' p)(p± we have 

1™ p 1.9? 2X 
M± = 2 F± " F± " T^^y 



F±xF± |</| 2 £pp + {Xn/a)&±af»l> 



2i 



X-Y 
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hence 



v±(x) = ^ p + ^ z ± ^ =► 4 = 1. 
Using ap = s/XY, we obtain (105). To derive (106), use the identities 




Note that u± is time-independent because the pulse factors \g'\ 2 in S± and u± 
cancel. Therefore the integral curves of v± are straight lines in M 3 propagating 
at the speed of light. Accompanied by their fields (E±, B±), these can now be 
justifiably interpreted as light rays. 

Whereas the physical interpretation of the null congruence associated with the 
scalar wavelet \P was vague because scalar (longitudinal) waves would seem to 
propagate in the direction orthogonal to their wavefronts, so their congruence 
could not be twisted, the identical null congruence defined by v± has a clear 
physical interpretation: the spinning disk T> emits light rays whose trajectories 
are given by 



x ± (£) = Po + tv° ± , Po eV, t>0 



u c V a2 - Pi ~ , 9 o t 
where v± — Sgn z — z ± — cp. 

a a 

The rays with p = form vertical vortex jets along the z-axis, while the 
rays emanating from any circle of radius < po < a generate the hyperboloid 
H 7)0 U H-^, with opposite rulings on the upper and lower halves. The rays 
coming from the circle po = a are tangent to that circle, forming the degenerate 
hyperboloid Ho. 

To find the complex congruence of F±, define the complex electromagnetic 
energy and momentum densities 

u= ^(E 2 + B 2 ), S = ExB. (107) 

Letting 

L = cos $ + k, M = A cos d + p 
and (90), we obtain 

1 + A 2 g 2 g 12 



(L 2 +M 2 ) (108) 

12 



2 C 4 P 2 

9 4(L 2 + M 2 )C +^f(L + XM)d + %f 
9 C P C P 
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Now assume that F± is null. Then A = and 



L 2 + M 2 = n 2 + p 2 = q+q-, L + AM = -q±, M — XL = ±iq±, 



so 



9 s 9 s 55 ~ 

u = q+q-'—r, S = q + q--—C-q±-^-^p±- 
p z p z Q^p 



Hence the complex energy propagation velocity is 



_ s s -\ 99 - 
5 C 2 



(109) 



(110) 



We shall suppress the helicity index ± and write v for notational simplicity. 
Thus 

v 2 = <f = l^c 2 , 
so the corresponding complex energy flow velocity 4-vector 

S" = (*,!), w M = (-v,l) 



defines a null congruence in C : 



= 0. 



To show that this congruence is twisted, consider the complex 1-form 

2 

5 P 

v = v u dx^ — dt — v • dec = dt — d£ + qZ ——r dtp± , 

+ 5 C 2 

where we have used (96). Since r = t — is and s is fixed, 
df — d£ = di~£ where T£ = r — £ 
is the complex retarded time. Letting 

and recalling that p/£ = sin$, we have 

S = dr^ + /i sin $ di? ± i/i sin 2 i? d</> 



The exterior derivative [F63] of u is given by 

dv = sin i? d/i A d?9 ± z sin 2 i? d/i A d</> ± ih sin(2tf ) di9 A dej), 



(111) 



(112) 



(113) 



(114) 
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hence 

v Adv = ±ih sin(2#) dr c A d$ A &<j>. 

Since g(r^) is analytic, its zeros are isolated and so are those of ft(r^). It follows 
that the congruence is twisting wherever sin(2z9) ^ 0. But 

sin(2i?) = 2sintfcostf = ^ = => p = 

C 

since z = z — ?a 7^ for all a; if a > 0. The complex congruence is therefore 
twisting everywhere outside the z-axis. Recall from (37) that the real congruence 
defined by E± ± iB± was twisting everywhere outside the xy-plane. 

Remark. The complex electromagnetic energy-momentum density is related 
to the real fields (E±, B±) as follows: 

u = l - (E + -E_+ B+ ■ B_) + 1 -{E_-B+-E + - B_) 

1 i 
S = -(E + x B +E x B + ) + -(E + x E_+B + x B_) 

They represent a Hermitian cross- correlation between the fields (E+,B+) and 
(E_,B_) whose physical significance, if any, remains to be investigated. 



11 Relation to Kerr-Newman black holes 

The Kerr-Newman solution of the Einstein-Maxwell equations in general rela- 
tivity can be described in a simple way as follows [Wik-KN]. 

• The spacetime metric ff~ v of the solution is given in terms of flat Minkowski 
spacetime coordinates (x,y,z,t) by 

9% = V^+fkpt (H5) 

where rj^ = diag(l, —1, —1, —1) is the Minkowski metric, the scalar function / 
is given by 

1 = c *™r-f) (11 6) 

where G is Newton's gravitational constant, and 

,4. ., , , , / rx =F ay ry ± ax z \ ,„ . 

K = ( fe ±- 1) ^ k ± = ( > r ) • (117) 



This describes a spinning black hole with mass M, charge Q, and angular mo- 
mentum L = mcaz. It can be checked that k 2 ± = 1, hence k^k^ = and k^ 
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defines a null congruence. Since k± is time-independent, the rays of this con- 
gruence are straight lines in Minkowski space. They become geodesies in the 
curved space defined by g^ v . 

• The electromagnetic field of the Kerr-Newman black hole is derived from 
the 4-vector potential 

Qr 3 ± 

where we have suppressed the helicity index ± in and F^. The parameter 
a above coincides with our parameter a, which represents the magnitude of the 
imaginary displacement of the point source. However, the parameter r does not 
signify the Euclidean distance |as| in M 3 . Instead, it identical to our parameter 
£, representing the real part of the complex distance C = £ — if]- Since az = £77 
by (20), / and simplify to 



Theorem 7 The 4-vector kjj; defining the null congruence of the Kerr-Newman 
solution (115) is identical with the 4-velocities (32) and (105) defining the 
null congruences of the scalar wavelet ^ (12) and the coherent electromagnetic 
wavelets F± (104), respectively. 



Proof: Since 



=Fayx ± axy — ±apq> and - = — , 



(117) gives 



£p±ap4> 7] £p ap ^ , ??- 

fe ± = ~C2~I 2~ + " Z = 7^- 2 P ± ~C2~, 2 <P + " Z - 



By (20) we have ap = \J + a 2 \J a? — n 2 , hence 




fc± = 

a 



However, note that the electromagnetic field F± of our coherent wavelets does 
not coincide with that of the Kerr-Newman solution. 14 In particular, F± is 



14 The complex representation of F M „ corresponding to F± is given in 4D spacetime notation 
by the field 

where *F I1V = \e^ va pF a ^ is the (Hodge) dual of F^ v . Since ** = —1, F M „ is self-dual: 

±i*F±=F±. 
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pulsed while F^ u in (118) is stationary. Both fields conform to the oblate 
spheroidal geometry of the Kerr-Newman spacetime. It remains to be seen 
how our coherent wavelets fit in with the theory of spinning black holes. 

Could a pulsed, radiating version g^ v of the Kerr-Newman metric (115) ex- 
ist which complements F± ? If so, this would give the following sequence of 
extensions: 

Scalar wavelet iff — > EM wavelet F± — > gravitational wavelet g^ v . 

Remark. Recall the general form of the axisymmetric complex vector potential 
A complementing iff : 

A(x,t) = iff(x,t)w(x), w = C + -(cos-& + k)S + -(Acostf + n)4>. 

P P 

The complex 4-potential (A, if>) for the general EM wavelet F± = E ± iB (not 
necessarily null) can thus be written in the form (119) as 

An(x,t) = <!>(x,t)wp(x) where w fJ , = (w,l). 

defines a complex null congruence if and only if w 2 = 1, or 

(cos $ + n) 2 + (A cos i9 + p) 2 = ^ 1 + A 2 = k + Xp = k 2 + = 0. 

The conditions for to define a complex null congruence are therefore 

~ c 

A = =pi and k = ±iw, hence w = £ H — (cost? + k)<p=c. (120) 

P 

By (91), this implies p± = q± = 0, hence F± = 0. 

That is, Wn defines a complex null congruence if and only if the Geld F± vanish 
identically. This can be checked directly by computing 

E = -^-C — — (cos$ + k)¥>t and B = =F«E. 
C 2 P 

This is an example of a 'pure gauge' field in the complex sense of (99), where 
A^O but F± = 0. 

Nevertheless, the null fields F± defined by A M with A = =F« and K ^ ±i/i define 
real and complex twisting null congruences (Section 10), and this the important 
fact since potentials have no direct physical significance. 
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